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Abstract. A study of the ultraviolet continuum variabil¬ 
ity (Paltani & Courvoisier 1994) has shown that the rel¬ 
ative variability of quasars and Seyfert galaxies decreases 
when the luminosity increases. The spectral information 
included in the lUE spectra allows us to study this de¬ 
pendence in the rest frame of the objects. The trend is 
strengthened by the general property that active galactic 
nuclei vary more at short wavelengths than at long wave¬ 
lengths in the ultraviolet domain. The scatter observed in 
all other studies is still present. An important part of this 
scatter may however be explained if one tries to estimate 
the uncertainties on the variability due to the sampling. 

We discuss the variability using the concept of discrete 
events. The trend between variability and luminosity is 
described by a power-law with an index —0.08, which is 
incompatible with the power-law of index —1/2 predicted 
by the most general discrete-event models in which the 
change in average luminosity is due to differences in aver¬ 
age event rates exclusively. Several biases are investigated, 
but we conclude that the —1/2 index is definitely inconsis¬ 
tent with the data. A flat relationship is however possible, 
if some bias has been underestimated. 

We propose different ways whereby discrete events may 
produce a different variability-luminosity relationship: by 
changing the luminosity or the life time of the events, or 
by introducing interdependence between the events. The 
latter possibility cannot produce a satisfactory relation¬ 
ship. Using the former possibilities, we do not find any 
“natural” explanation for the variability-luminosity rela¬ 
tionship in the context of discrete-event models. This is 
possibly an indication that explanations in which variabil¬ 
ity is not expressed in terms of discrete events should be 
favoured. 
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1. Introduction 

The statistical properties of the ultraviolet variability 
of active galactic nuclei (AGN) have been investigated 
by Paltani & Gourvoisier (1994) (hereafter PG94). This 
study, which was based on spectra obtained by the Inter¬ 
national Ultraviolet Explorer (lUE) from the ULDA ( Uni¬ 
form Low Dispersion Archive) database, has been done 
in the observer’s frame. In addition to the variability at 
2 000 A, we quantified also the change of variability with 
the wavelength, which showed that all classes of AGN have 
a larger variability at short wavelength than at long wave¬ 
length in the lUE domain. In this paper we shall consider 
only the objects that have been called Seyfert-like objects 
in PG94, i.e. Seyfert 1 galaxies, radio-quiet quasars, and 
low-polarization radio-loud quasars. 

The relationship between variability and luminosity 
in quasars and Seyfert galaxies has been investigated by 
many authors. Most of them found an anti-correlation be¬ 
tween these two quantities (Pica & Smith 1983; Gristiani 
et al. 1990; PG94; Hook et al. 1994), while others found 
a correlation or no correlation at all (Bonoli et al. 1979; 
Trevese et al. 1989; Giallongo et al. 1991). Rediscussion of 
these last three papers by Hook et al. (1994) seems to show 
that all the samples are compatible with the existence of 
an anti-correlation. 

All these studies have used estimates of the variability 
in the observer’s frame, as they can be much more easily 
obtained than those in the objects’ rest frame. Therefore 
the correct relationship between variability and luminosity 
is still unknown. This relationship is very important, as it 
can be compared with the predictions of some AGN mod¬ 
els, in particular the starburst model of Terlevich (1992). 
In this paper we establish this relationship in the objects’ 
rest frame, and we discuss some biases that can affect it. 
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We then check the compatibility between our result and 
the predictions of these models. 

We formalize in the most general way the concept of 
“discrete-event models”, i.e. models where the total light 
curve is the sum of a series of light curves, each one asso¬ 
ciated with one “event”, whatever its physical nature, and 
we compare their predictions with the above constraint. 
We give the general relationship between variability and 
luminosity for these models. We study two different ways 
by which this relationship can be modified. In the first 
one we assume that some of the event parameters may 
actually depend on the luminosity. In the second one we 
examine the effects of the introduction of interdependence 
between the events. 


2. Rest-frame variability in the ultraviolet 

2.1. Determination of the rest-frame variability 

In this paper, as in PC94, we use the lUE spectra from the 
ULDA version 4.0 database. This database contains nearly 
every low-dispersion spectrum obtained by the lUE satel¬ 
lite before 1992. The period covered by the ULDA archive 
is therefore about 14 years. The wavelength range covered 
by these spectra, combining short- and long-wavelength 
cameras, is from 1200 A to 3200 A. The large number of 
repeated observations and the spectral coverage of the lUE 
satellite make this satellite ideal for the determination of 
the AGN variability properties. 

In PC94 we estimated for each object the wavelength- 
dependent variability, cr/(A). In all this paper, as in PC94, 
“variability” should be understood as the standard devia¬ 
tion of the flux divided by the mean flux. For each object 
we approximated cr/(A) by a linear relationship: 


cr/(A) ~ Vx{af) ■ 


A - 2000 A 
1000 A 


+ O’/,2000 


( 1 ) 


O ’/,2000 is expressed in percents of the mean flux and 
Va(o’/), the slope of o’/(A), is expressed in percents of 
the mean flux per 1 000 A (cf. Fig. 1 of PC94). The values 
of O’/,2000 and Va(o’/) can be found in the Appendix of 
PC94. The definition of these parameters has only a de¬ 
scriptive purpose and does not have any predictive value 
outside the wavelength range studied. Using these two pa¬ 
rameters, we can derive the AGN variability at a given 
wavelength, 1 250 A in our case, in the rest frame of the 
objects: 


O'/;i250 - O’/,2000 + Va(0’/) 


1250- (1 -f z)- 2 000 
1000 


( 2 ) 


We feel that it is important not to extrapolate the vari¬ 
ability outside the lUE wavelength range. We therefore 
consider in this paper only the objects that have a red- 
shift smaller than 1.3. 


2.2. Uncertainties on the rest-frame variability 

It is important to estimate the uncertainty on o’j®® 250 ’ 
we shall be concerned in this paper with the relationship 
between variability and another parameter, namely the 
luminosity. The uncertainty on crj ®®250 origins. 

First it is due to the propagation of the uncertainties on z, 
o’/, 2000 j and VA(cr/) from PC94. As explained in PG94, the 
two last uncertainties are the uncertainties on the least- 
squares linear regression and do not take into account the 
effect of the sampling. A second term must therefore be in¬ 
cluded to allow for the effect of the sampling. We propose 
here a method that is an improvement of the discussion 
of Sect. 3.2 in PC94. 

The first step is the simulation of a light curve that has 
variability properties comparable to those of our objects. 
We assume that the AGN light curves are characterized by 
a broad Fourier power spectrum which follows a power-law 
with an index between —1 (flicker noise) and —2 (random 
walk). These has been verified for at least two Seyfert-like 
objects: 3C 273 (Kunkel 1972) and NGC 4151 (Longo et al. 
1996) (see also Press 1978). An example of a light curve 
simulated by a “random walk” can be seen in Fig. 1 of 
PC94. We rescale these simulated light curves to obtain a 
variability comparable to the ones observed in our objects. 
Then we draw n epochs from an uniform distribution on 
the period covered by ULDA and compare the variability 
estimated with these n points with the variability of the 
simulated light curve. Repeating this last step many times, 
we can estimate the uncertainty produced by the sampling 
from the dispersion of the estimates of the variability. 

As the indices are smaller or equal to -1, a low- 
frequency cut-off has to be chosen if one wants to avoid 
inhnite variability. On the basis of archival data of 3C 
273 (Angione & Smith 1985), it appears that variability 
on time scales longer than 20 years is not dominant. If 
this is valid for all AGN, a 14-year period should contain 
almost all the power of the Fourier spectrum. The effect 
of variations on longer time scales and of time dilation 
in high-redshift objects are discussed in Sect. 3.1 and 3.2 
respectively. 

Taking 3 values for the spectral index of the power 
spectrum (—1, —1.5, and —2), 3 values for the low- 
frequency cut-off ((5 yr)“^, (10 yr)“^, and (20 yr)“^), and 
3 values for the relative variability (0.05, 0.1 and 0.5), we 
parametrize the relative uncertainty on a, as a func¬ 
tion of n: ^(n) = CT 2 o • (n/20)® (where CT 2 o is the uncer¬ 
tainty for n = 20), which produces good approximations. 
We normalize — at n = 20 because this value is central 

a 

and is little affected by small deviations of s. It appears 
from the simulations that none of the 3 parameters modi¬ 
fies significantly the value of s, and s = —0.56 reproduces 
quite well the decrease of at least for n < 50. On the 
other hand, (120 depends on the 3 parameters, but the dif¬ 
ferences are only a few percents. 1720 increases when the 
index increases, the low-frequency cut-off decreases, or the 
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Fig. 1. crJ®j25o vs luminosity for Seyfert-like objects. The solid 
line is the power-law obtained with the BCES method, and 
has a slope of -0.078. The dashed line is the best fif with a 
power-law of fixed index —1/2. 

variability increases. We therefore use the following empir¬ 
ical relation for the uncertainty on the relative variability 
due to the sampling in the worst case, i.e. low-frequency 
cut-off=(20 yr)“^, and index= —1 (flicker noise): 

( 3 ) 

The best case makes (T 20 become (0.14 -|- 2 • ct^). In prac¬ 
tice, as the number of spectra is usually different in the 
short- and long-wavelength ranges, n will be the larger 
of the two numbers. For the few cases where the number 
of spectra is larger than 50, we use = 0.05. It is im¬ 
portant to note that the estimates of the variability follow 
rather well a Gaussian distribution as soon as n is larger or 
equal to 4 (the distribution is already marginally accept¬ 
ably Gaussian for n = 3). Thus we can estimate the total 
uncertainty on simply adding the variances of 

the two sources of dispersion. 

2.3. Correlation between variability and luminosity 

Figure 1 shows the variability ^ function of the 

luminosity for all the objects classified as Seyfert-like ob¬ 
jects in PC94 (the objects and their properties are given 
in the appendix of PC94). The Spearman’s correlation co¬ 
efficient is equal to -0.45. The probability to get such a 
high correlation coefficient for 68 couples drawn from two 
uncorrelated parent populations is about 10“^. We note 
that Spearman’s test is valid even if the variables have 


uncertainties or intrinsic scatter. A real correlation may 
be destroyed by the dispersion, but the probability that 
the dispersion generates a spurious correlation is very low. 
Therefore, with better measurements, the correlation co¬ 
efficient should in principle increase (in absolute value). It 
appears therefore unquestionable that a trend exists, and 
this trend is in the sense that the variability decreases 
when the luminosity increases. 

To quantify this trend, one must choose an analytical 
formulation. Because the very large dynamic of the x axis 
(~ 10^), the most natural one is a power-law relation in 
the form: 

'^/4250 ~^1250A (4) 

where L ° is the average luminosity density of the 

r0SL ] 

object at 1250 A in the rest-frame of the object. The in¬ 
dex Tj of the trend can in principle be estimated. However, 
because of the important uncertainties, one must be care¬ 
ful about the selection of the method. As we determined 
both X- and y-axis uncertainties, we used a method that 
takes into account the uncertainties on both coordinates 
(Press et al. 1992). This method is valid only for linear 
(affine) relationships; therefore we fit the logarithm of the 
values. Unfortunately, we have to estimate the variances of 
the logarithms of the distributions, which means that the 
uncertainties are not Gaussian anymore. We nevertheless 
assume that they follow closely enough Gaussian distribu¬ 
tions in order to apply the test. We obtain rj = 0.084. 
The probability that a linear relationship really holds is 
of the order of 10“^ (x^ — 104 with 66 degrees of free¬ 
dom). As our uncertainties are probably overestimated, 
some intrinsic scatter around this relationship does exist, 
but a large part of the dispersion can be accounted for by 
the uncertainties. The linear function with a slope fixed 
to —1/2 that minimize the x^ has been plotted on Fig. 1. 
It is clearly unadapted to the data (its x^ is larger than 
1000). The value of the slope obtained using the usual 
least-squares method (which neglects the uncertainties on 
both axes) is —0.078, which is very close. 

The scatter makes the uncertainty on the slope of the 
relationship meaningless. Fortunately a new method has 
appeared that is not only applicable when uncertainties 
on both axes exist, but also when there is some intrinsic 
scatter. This method, called BGES (Akritas & Bershady 
1996), gives a slope r] = —0.078 ± 0.16. The agreement 
between all these methods reinforces the confidence in the 
value, and we use in the following: 

r] = -0.08 (5) 

The BGES method is, as the authors claim, very con¬ 
servative. Although ?7 = 0 is formally possible, Spearman’s 
test, which is much more efficient in detecting the exis¬ 
tence of a slope, clearly rules out this hypothesis. 

It is interesting to note that Hook et al. (1994) have 
found indices of —0.06 (their model A) and —0.13 (their 
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model B) for the same relationship, which are both com¬ 
patible with our result (but see Sect. 3.3). Their quasar 
sample extends far beyond our sample in redshift and 
in luminosity and does not comprise any Seyfert galaxy. 
Moreover, as we discuss in the next section, they do not 
allow for several biases. 

3. Biases in the variability—luminosity relationship 

We have seen in Sect. 2.3 that the slope of the variability- 
luminosity relationship, although statistically different 
from 0, is small. Therefore, we can ask whether this trend 
could be due to a bias. On the other hand, it is also im¬ 
portant to ask whether the slope of the relationship might 
be lessened by systematic effects. 

3.1. Very-long-term variability 

The ULDA 4.0 data span about 14 years. It means that 
variability on time scales longer than ~ 20 years will not 
contribute (or only partly) to our estimation of <7/(A). If 
only low-luminosity objects vary on very long time scales, 
the observed (j{L) relationship would be flatter than the 
intrinsic one. On the other hand, the decrease of variability 
with luminosity may be introduced by bias if only high- 
luminosity objects vary on very long time scales. 

Very-long-term variability can have two manifesta¬ 
tions: either slow, steady variations or short, intense flares 
occuring with a rate of around 1 per century, or less. In 
any case, this very-long-term variability implies that the 
Fourier power spectrum of the AGN light curves should 
contain a non-negligible amount of power at frequencies 
smaller than (20yr)“^. We can estimate this power for 
two limiting cases assuming a given intrinsic cr(L) rela¬ 
tionship: 

— The variability is independent of the luminosity. In this 
case the more luminous the object, the larger the unob¬ 
served proportion of power in the Fourier spectrum at 
very small frequencies. Assuming that only variability 
on short time scales exists in low-luminosity objects, 
it implies that short-time-scale variability represents 
only 40 % in high-luminosity objects. The variabil¬ 
ity of 3C 273 from optical archival data since 1887 
(Angione & Smith 1985) is about 16 %. Its variability 
on a lO-year basis is 12 % (work in progress). In the 
Fourier spectra presented by Angione & Smith (1985), 
much power seems to be contained between (20 yr)“^ 
and (10 yr)“^, which could easily explain the differ¬ 
ence, but the increase in variability is much smaller 
than the required factor 2.5. The power contained be¬ 
tween (100 yr)“^ and (20 yr)“^ seems negligible. Even 
though it is a priori possible that still much power ex¬ 
ists at frequencies below (100 yr)“^, the small increase 
of power between (100 yr)“^ and (20 yr)“^, together 
with the necessity of the disappearance of this power 


in low-luminosity objects, makes this possibility rather 
improbable. 

— The intrinsic a{L) relationship has a —1/2 slope, which 
is the slope predicted by many discrete-event mod¬ 
els (see Sect. 4). In this case the variability of low- 
luminosity objects must have been underestimated by 
a factor about 30. Assuming that the variability of 
high-luminosity objects is correct, which means that 
only variability on short time scales exists in these ob¬ 
jects, the luminosity of low-luminosity objects should 
be dominated in average by very-long-term variations. 
As short-term variability is not much smaller than 100 
% in most low-luminosity objects, it would imply that 
all these objects are in phases where the very-long-term 
variations have a very small luminosity. This is statis¬ 
tically impossible if the very-long-term variations are 
slow and steady. Short intense flares can also be dis¬ 
missed, as a sufficient variability can be attained only 
if the average luminosity is dominated by the very rare 
flares by a large factor. Such flares have never been ob¬ 
served. Taking into account the intense monitoring of 
several Seyfert 1 galaxies, the hypothetical flaring state 
can be active only for a very small fraction of the time. 
Consequently, the flare intensity should be extremely 
large, should have been detected in the same way as 
supernovae are detected. 

In conclusion, although we cannot completely exclude 
any important effect on the variability-luminosity rela¬ 
tionship, it appears more probable that such biases are 
marginal. But if the variability-luminosity relationship is 
affected, it is certainly in the sense that the trend is arti¬ 
ficially produced rather than in the sense that the —1/2 
relationship is flattened. 

3.2. Time-dilation effect 

If light curves of AGN are stochastic time series that have 
Fourier power spectra well described by a power-law with 
indices between —1 and —2, the total variability of the ob¬ 
jects increases with the length of the sampling until this 
length becomes comparable to the inverse of the frequency 
of the low-frequency cut-off. Because of the time dilation, 
high-redshift objects are potentially observed for a smaller 
time than low-redshift ones. Let us assume that the “stan¬ 
dard” index is —2, which corresponds to a “random walk”. 
If x{t) is such a time series, its structure function Sx{t) 
(Rutman 1978; Simonetti et al. 1985) is linear: Sx{t) ~ r. 
Larger values of the power-spectrum index make Sx(t) 
increase more slowly. 

The value of the structure function is also related to 
the variance of the light curve: 

S'2,(to) ~ 2 • Var(xrJ (6) 

where Var(xTo) is the variance of a portion of the time se¬ 
ries x, whose length is tq. However, because of the redshift 
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of the sources, one does not measure the correct length 
of the time series tq, but the observer’s length, which is 
r = To • (1 -b z), where z is the redshift. Therefore, assum¬ 
ing that the observer’s lengthes of the time series are all 
identical, one has to modify the observed variability: 

2Var(a;^) ~ Sx(t) = (1 + z)Sx(to) ~ 2 (1-b z) Var(xro)(7) 

which gives: 

— ^X,To ' “b "z (8) 

This correction is an upper-limit, because it assumes 
that the power spectrum extends down to frequencies 
much smaller than (14 yr)“^. If one applies this correc¬ 
tion to our data, the anticorrelation still holds with a con¬ 
fidence of 99.8 %, and the slope obtained by the BCES 
method is —0.058 ± 0.15. Therefore, it seems that the de¬ 
crease of the variability with the luminosity cannot be 
entirely explained by this effect. 

3.3. Existence of a correlation between variability and red- 

shift 

Luminosity and redshift are always very well correlated in 
flux-limited samples, since only the most luminous objects 
can be observed at large redshift. This leads to a corre¬ 
lation between variability and redshift. The question that 
arises is to determine whether this correlation is entirely 
due to the bias, or whether there also exists an intrin¬ 
sic correlation. Hook et al. (1994), using subsamples with 
small luminosity range, have concluded that the first al¬ 
ternative is the correct one. 

Two effects can produce a correlation between variabil¬ 
ity and redshift. The first is evolution; unfortunately, we 
are much too far from a clear picture of how AGN work to 
guess anything about its effect, if any, on the variability- 
redshift correlation. The second effect is much clearer; in 
PC94, we showed that the variability increases with the 
frequency. As the observations are usually performed at a 
wavelength fixed in the observer’s frame, this should nat¬ 
urally lead to a positive correlation between variability 
and redshift. Di Clemente et al. (1996) found that this ef¬ 
fect alone can account for the positive correlation between 
variability and redshift, which means that there is no need 
to invoke evolution effects. 

To allow for this effect, Cid Fernandes et al. (1996) 
have fitted the data of Hook et al. (1994) with a function 
of two variables, the luminosity and the redshift, and three 
parameters, a, 5, and c: 

logo-= a • logLs-b 5 • log(l-b z)-b c, (9) 

where Lb is the mean luminosity in the optical B band 
and z is the redshift. The parameter a is is equivalent 
to our parameter ry; b takes into account the increase of 
variability with the frequency; and c is a normalization 
parameter. Although evolution effects could in principle 


affect parameter b, this is not considered by Cid Fernandes 
et al. (1996), because their preferred model makes use of 
independent events, identical for all objects (see Sect. 4). 

The results of their study are: —0.6 < a < —0.2 and 
0.1 < b < 1.0 (in one particular case; however other cases 
are qualitatively identical). It is in clear contradiction with 
the results cited above, in particular ours. The explanation 
of the difference between this result and the one of, for 
instance. Hook et al. (1994) lies essentially in two points: 

1. Cid Fernandes et al. (1996) take into account the fact 
that AGN vary more at shorter wavelengths. Hook et 
al. (1994) neglect entirely this effect after checking that 
it is not important in their sample. 

2. Cid Fernandes et al. (1996) estimate the uncertainties 
on the luminosity, and take them into account in the 
fits. La Franca et al. (1995) have shown that the slopes 
are biased towards a smaller value when the errors on 
the X axis are discarded. 

Since our result is very close from the one of Hook et al. 
(1994), who neglect these two biases, further discussion is 
required. 

First, we note that the increase of the b parameter 
has an effect very close to the decrease of a. This implies 
that solutions with a ~ —0.5 are indeed allowed by Cid 
Fernandes et al. (1996), but it requires large values of b 
(~ 1). Although such values may be allowed in some cases, 
it cannot be guaranteed that most objects follow this be¬ 
haviour. The determination of the rest-frame variability 
removes the difficulty encountered by Cid Fernandes et 
al. (1996) about the value of b. In our case, no hypothesis, 
no extrapolation is needed. As evolution is probably not 
a problem (the more so as all our objects have redshift 
smaller than 1.3), no correlation between variability and 
redshift is expected, apart, perhaps, from the time-dilation 
effect discussed in Sect. 3.2. 

We can object to the second point that we also took 
into account the uncertainties on the x axis to obtain the 
slope. On the other hand Cid Fernandes et al. (1996) did 
not estimate the uncertainties on the variability, which, 
in our case, are much larger than those on the luminos¬ 
ity. They obtain steep values for the variability-luminosity 
relationship only in the cases where the uncertainty on 
the luminosity is very large (which is, by the way, an ex¬ 
cellent illustration of the bias described in La Franca et 
al. (1995)). Moreover, only the BCES method takes cor¬ 
rectly into account the existence of intrinsic scatter in the 
data. Once again, we took advantage of the spectral ca¬ 
pabilities of lUE: while Cid Fernandes et al. (1996) have 
to derive /c-correction based on unknown spectral indices, 
our luminosities are automatically calculated at 1250 A 
in the rest frame of the objects, again without any hy¬ 
pothesis or extrapolation. We note that the —1/2 slope 
is compatible with our data at the 1 a level if the un¬ 
certainties on the luminosity are multiplied by a factor 
18 {rj = —0.14 ± 0.36), which is clearly an overestimate 
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(moreover, the BCES method has very conservative un¬ 
certainties). 

As a consequence we infer that the difference in the 
result is due mostly to the assumptions that they had to 
make to correct from the wavelength dependence of the 
variability and to obtain rest-frame luminosities. 


analytical development is made in Appendix A and the 
result is given in Eq. (A5). If one makes the assumption 
that the event rate is independent of the parameters defin¬ 
ing the events, one finds that (from Eq. (A6)): 

Var(L)~7V (13) 


4. Variability from discrete events 

4.I. The concept of the discrete-event models 

It is very convenient to think of the variability as due to 
the superimposition of perturbations on a constant source. 
The luminosity density at frequency n can then be written: 

L^{t) = C„e^^i{t - U), ( 10 ) 


where is the constant state of the object and is the 
light curve produced by one event (the index i indicates 
that the event light curve may be different for each event), 
and ti is the “birth time” of the event. It is a generaliza¬ 
tion of a shot-noise process. This representation can apply 
to the blue-bump of the Seyfert-like objects, as the vari¬ 
ations are quasi-simultaneous from the ultraviolet to the 
optical (Courvoisier & Clavel 1991). Moreover it is very 
general, and can be a reasonable approximation of many 
models, from accretion disk models, e.g. the “blob” model 
(Haardt et al. 1994) or shocks in the disk (Chakrabarti & 
Wiita 1993), to supernova explosions (Terlevich 1992) and 
relativistic star collisions (Courvoisier et al. 1996). The 
interesting point of the discrete-event models is that the 
variability properties of a set of object is subject to strong 
constraints as soon as one puts limitations on Eq. (10). 

To demonstrate analytically these constraints, we de¬ 
velop the following formalism: The mean luminosity den¬ 
sity L^{t), which can be written (in the following we omit 
the ‘V” index, the terms “luminosity” and “energy” must 
be understood as “luminosity density ” and “energy den¬ 
sity”, and we assume that the studied systems are com¬ 
pletely stationnary): 

1= lim — / L{t) dt = N -E + C, (11) 

where N is the event mean rate, Ei = ei{t — ti) dt is 
the total energy released by the event (at frequency v) 
and E is the mean of Ei over all Ts. Eq. (11) is valid even 
if the ei{t) differ from each other. 

It is much more complicated to obtain the variance of 
L{t). Its exact definition is: 

Var(L) = lim E f {L{t)-L)^dt (12) 

1 —^oo ZJ. j _ 'j' 

Using the hypothesis that the light curve is given by 
Eq. (10), it is possible to transform this equation. The 


The relative variability, i.e. the ratio between the standard 
deviation (the square root of the variance) and the mean, 
is therefore: 


cr(L) 


Vn 

N-E + C 


(14) 


The strongest limiting case on Eq. (10) has already been 
discussed by Pica & Smith (1983): the whole luminosity is 
due to the events {C = 0) and all events are independent 
and identical. Eq. (14) shows that (t(T) is proportional to 

N , i.e. proportional to L (by Eq. (11)). However, 
Eq. (14) is much more general; in particular the events 
may be completely arbitrary, and the —1/2 index will still 
be preserved. C = 0 is not even necessary. The —1/2 in¬ 
dex is also obtained if C is a fixed fraction of the total 
luminosity {C = x ■ L, 0<a;<l). Moreover, if C is con¬ 
stant for all objects (e.g. if C is the contribution of the 
underlying galaxy), we still find the —1/2 index as soon 
as C is negligible compared to the events’ luminosity. In 
conclusion, we see that the only requirements to preserve 
the index —1/2 is that the parameters of the events are 
drawn from the same distributions in all objects (i.e. the 
distributions of the parameters are independent of A^). 

Pica & Smith (1983) concluded therefore that such 
models can be rejected independently of the nature of 
the events (although they did not prove this rigourously). 
They have found that the —1/2 index could be dismissed 
with a confidence larger than 99 %. PC94 and Hook et 
al. (1994) have reached a similar conclusion. However in 
all these studies the variability have been made in the 
observer’s frame. In this paper, the variability used to de¬ 
termine the relationship with the luminosity is corrected 
for the redshift effect, and here again the value —1/2 can 
be rejected with a large confidence. 

Terlevich (1992) proposed that the AGN could be pow¬ 
ered by independent supernova explosions. In this model, 
the events (the supernovae) are identical and independent, 
and C' = a;-L, 0<a;< lis satisfied, because a fixed 
fraction of the total luminosity is emitted by normal (con¬ 
stant) stars {x is larger than 0, but depends strongly, of 
course, on the considered wavelength). As we have seen, 
the index of the relationship between luminosity and vari¬ 
ability remains equal to —1/2. 

Even if a general trend is predicted, we see that much 
room is left for the observed intrinsic scatter in the terms 
inside the integrals in Eq. (A6). Eor instance the distri¬ 
butions of the luminosities of the events may differ from 
object to object without being necessarily related to the 
total luminosity of the object. 



S.Paltani & T.J.-L.Courvoisier: Rest-frame variability of AGN and discrete-event models 


7 


As the observations clearly prohibit the so universal 
— 1/2 index, we could ask whether the release of one of 
the two properties of the events (“identity of the distribu¬ 
tions of the event parameters” and “independence”) might 
produce the observed behaviour. 


J^.2. Non-identity of the distributions of the event param¬ 
eters 

The variability-luminosity relationship can be modihed to 
match any given relationship if one adjusts in an ad-hoc 
way the distributions of the event parameters. We examine 
only two very simple cases that can be intuitively expected 
in many discrete-event models. 


4.2.1. Scaling of the event luminosity 

A very simple modihcation that breaks the “identity of 
the distributions of the event parameters” assumption is 
that the energy released by the events scales with the total 
luminosity of the object, while the events are still indepen¬ 
dent from each other. We can estimate the effect of such a 
scaling on the variability-luminosity relationship using the 
notation of Eq. (Bl) and the results of Eqs (B4) and (B3) 
(see Appendix B) in the case C = x-L, 0<x<l. The 
scaling is written E = k-Eo ^ L°‘, which gives N ^ L ^. 
Therefore Eq. (14) becomes: 


a{L) 


I 


-(a-l)/2 


(15) 


If no correlation exists between variability and luminosity 
(?7 = 0), then we must have: a = 1. If the rj = —0.08 index 
is correct, one obtains: 


4.2.2. Scaling of the event stretching 

We call “life time” of an event any time that is mathemati¬ 
cally associated with the event. For instance, the life time 
of a “box” event is the duration of the non-zero state; 
for an semi-exponential event, the life time may be the 
e-folding time. We may think that the event life time in¬ 
creases with the total luminosity. This may appear natu¬ 
rally in many models. For instance, in case of reprocessing 
on an accretion disk, the size of the reprocessing area (con¬ 
tributing at a given frequency) very certainly increases 
with the luminosity of the X-ray source(s). We call the 
possibility to modify the event life time the “stretching” 
of the events. 

We introduce in the mathematical formulation of the 
events a stretching factor £ (corresponding to “p 2 ” in 
Eq. (Bl)). We assume that (. scales with the total lu¬ 
minosity of the object, i.e. £ = 6 ■ ^ (see Ap¬ 

pendix B). From Eqs (B4) and (B3) we obtain (valid if 
C = X ■ L, 0 < a; < 1): 

„(L)_ ,17) 

L 

If ?7 = 0, we must have /3 = — 1; with our value for rj we 
find: 

(3 ~ -0.84 (18) 

These values are not satisfactory, because one would 
rather expect that the stretching factor increases with 
the luminosity, and not decreases. It appears therefore 
that this mechanism (alone, at least) cannot explain the 
variability-luminosity relationship. 

4.2.3. Scaling of both the event luminosity and the event 
stretching 


a-0.84 (16) 

This scaling cannot be applied to models based on super¬ 
novae, as the energy released by a supernova is essentially 
fixed. On the other hand, it is reasonable to assume that 
the luminosity of an event - of unspecihed nature - re¬ 
lated to the existence of a black hole scales - in a way 
or in another - with its mass. For instance, Courvoisier 
et al. (1996) have proposed a model where such a scaling 
with luminosity could be rather natural. In this model, the 
events are collisions between stars in a dense stellar clus¬ 
ter surrounding a supermassive black hole, which liberates 
the stars’ relativistic kinetic energy. For a given cluster, 
an increase of the black-hole mass will increase the energy 
liberated by the events. 

Generally speaking, however, the value of a shows that 
the event luminosity is not trivially related to the total 
luminosity. This may be a major difficulty for the discrete- 
event models. 


We can also combine the scaling of the event luminosity 
and the one of the event stretching. Again, the justihcation 
of these scalings is that they seem natural to us. Eqs (B4) 
and (B3) give here: 


<j{L) 



l/2+(a-/3)/2 


(19) 


If we assume that (3 = 1, which means that the stretching 
factor of the events is proportional to the total luminosity, 
our result requires: 


a - 1.84 (20) 

or, if a = 1, i.e. the event luminosity is proportional to 
the total luminosity: 

/3-1.84 (21) 


If we impose a > 0 and /3 > 0, which is physically the 
most appropriate, we Hnd that a must be > 0.84, because 
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j3 ~ a —0.84. If ?7 = 0, we must have a > 1, and /3 = a—1. 
The models are at the moment not sufficiently developed 
to predict any value for a and j3. We can however remark 
that the constraints on a and (3 are in no case intuitively 
“natural”. 

^.3. Interdependence of the events 

4.3.1. Pseudo-Poisson process 

It has often been argued that a different relationship be¬ 
tween variability and luminosity could be obtained if the 
events were not independent from each other. For instance, 
one can imagine a cluster of explosive charges. The charges 
may explode spontaneously (producing an “event”), or the 
explosion may be triggered by the close explosion of an¬ 
other charge. It is easy to imagine a mathematically anal¬ 
ogous situation in AGN. 

In this picture, each spontaneous event (not triggered 
by another event) is followed by a sequence of “child” 
events. Therefore one can replace the spontaneous event 
and its “child” events by one “super-event”. The normal¬ 
ization and shape of this “super-event” are not constant, 
but are drawn from complex distributions. 

It appears therefore from the discussion in Sect. 4.1 
that the —1/2 index cannot be avoided by any process 
following this picture. For instance, if all the events are 
independent and if the “parent” event has a probability 
(f to produce a simultaneous “child” event, which in turn 
has a probability tp to produce another “child” event, and 
so on, the variability at a given total luminosity (if (7 = 0) 
is multiplied by: ^/1 + — (p. 

4.3.2. Non-Poisson process 

The reason why the —1/2 index always holds in the pro¬ 
cesses dicussed above is that the process remains Poisso- 
nian in essence. However no clear physical picture of non- 
Poisson process can be proposed at the moment. We can 
nevertheless propose a mathematical non-Poisson process: 

In a system where events occur according to the Pois¬ 
son law, the probability Fo(G) that no event occurs be¬ 
tween t = 0 and t = ti follows the equation: 

^ • ^0, (22) 

with the boundary condition Poit = 0) = 1, p being the 
mean number of events between t = 0 and t = 1. The 
solution of Eq. (22) is simply: 

Po(f)=e-^‘ (23) 

We can modify Eq. (22), so that the mean time with¬ 
out event decreases when new events occur; this reflects 
the interdependence of the events. We write: 

M + £» • ^ fit-tkU ■ Po, (24) 

fc=i / 


where g is the “response” to the occurence of an event - 
whatever its origin n is the number of events that have 
already happened, and f{t) is the time distribution of the 
influence of an event, with f{t) dt = 1. g should be 
smaller than 1, as a larger value would produce a diverging 
“chain reaction”. 

It is possible to solve Eq. (24) numerically for each new 
event and to draw lists of successive events for a given 
choice of p and g. If we choose an event light curve, we 
can construct the total light curve, and deduce the lumi¬ 
nosity and the variability of the virtual object. In all the 
simulations, the distribution function f(t) has been taken 
to be uniform between 0 and I (arbitrary) unit of time 
(A) and the event light curve was a simple box function. 
For a set of g chosen between 0 and 0.96, we calculated 
different light curves, making p, the rate of independent 
events, vary from 10“^ to as many events per A as pos¬ 
sible (the limitation is due to the calculation time, which 
increases with the rate of events). Examples of light curves 
are shown in Fig. 2. 

Figure 3 shows the variability-luminosity relationship 
obtained from the simulations for different values of g. In 
all curves the behaviour is comparable; one can distinguish 
3 regimes: 

— The sub-critical regime occurs when the mean rate of 
events is small. The variability-luminosity relationship 
has the slope predicted by pure Poisson statistics, in¬ 
dependently of g. The events are too rare to affect sig¬ 
nificantly Eq. (24). The light curves (Fig. 2a) do not 
show clear structures (as is expected from pure Pois¬ 
son process). 

— The critical regime occurs when the mean rate of 
events is around I per A (slightly dependent on g). 
The variability decreases less rapidly, or even increases 
for large values of g. The larger g, the larger the de¬ 
viation from the —1/2 index. In this regime, the event 
density is right in the range that produces “chain reac¬ 
tion” . The light curves show very characteristic, flare- 
like structures due to the interdependence of the events 
(e.g. around t = 3100 A in Fig. 2b). 

— The super-critical regime occurs when the mean rate 
of events increases further. The variability-luminosity 
relationship has again the slope predicted by pure Pois¬ 
son statistics, but with a larger normalization than in 
the sub-critical regime. It seems that the “chain reac¬ 
tion” is saturated by the large event density, and the 
light curves (Fig. 2c) have essentially lost the char¬ 
acteristic structures, probably because of the “over¬ 
crowding” of the events. The increase of variability in 
this regime as a function of p is compatible with the 
multiplication by a factor ^/1 -f g/— g. 

Even though this model can produce a variability- 
luminosity relationship different from the Poisson case, 
none of the curves in Fig. 3 can be adjusted to the 
data in Fig. 1. Moreover, the physical feasability of this 
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Fig. 2. a-c. Portion of light curves simulated with g = 0.9 in 
the three different regimes: a sub-critical, L = 0.1 events/A, b 
critical, L — 2.2 events/A, c super-critical, L = 72 events/A. 
d Light curve simulated with £< = 0, and the same mean lumi¬ 
nosity as in panel b 



Log Luminosity [ Events/A J 


Fig. 3. Variability as a function of the luminosity (expressed 
in number of events per A) in the simulation of interdependent 
events for different values of g. The £< = 0 curve has a —1/2 
slope. 


model cannot be asserted. Other mathematical formula¬ 
tions might yield more appropriate results. However, the 
physical feasability might be a major difficulty in most 
(if not all) of them, because (pseudo-)Poisson processes 
appear much more natural. 

5. Conclusion 

The trend between the variability and the luminosity seen 
by many authors is confirmed when one estimates the rest- 
frame variability. The most probable index of the relation¬ 
ship is —0.08. It is compatible with results obtained using 
high-luminosity quasars. This is an additional argument 
for the continuity between the Seyfert galaxies and the 
most powerful quasars. 

We have investigated the effects of several biases, and 
we consider our result as very probably correct. An in¬ 
trinsic ry = —1/2 index, which is what is generally ex¬ 
pected from discrete-event models, would be possible only 
if there exists very large and yet undiscovered biases; but 
this seems highly unlikely. The other extreme, which is 
that the variability does not depend on the luminosity, is 
more plausible. However the difficulty in producing such 
a relationship with discrete events is not significantly al¬ 
leviated. 

The easiest way to produce a relationship with an in¬ 
dex different from —1/2 is to assume that the average 
event luminosity increases with luminosity. In this case, 
it would mean that the event luminosity must increase 
only about 15 % more slowly than the total luminosity. 
The scaling of the event stretching factor may also con¬ 
tribute to the relationship. The relationships given here 
may provide useful constraints for any kind of discrete- 
event model. 

The possibility that the birth time of the events does 
not follow a Poisson distribution has also been investigated 
and cannot be completely excluded. However the system 
that we have explored cannot (alone) produce the cor¬ 
rect relationship. Anyway a non-Poisson system appears 
somehow unnatural, and the physical situation at its ori¬ 
gin might be difficult to find. 

The predictions of physical models that can be ex¬ 
pressed in terms of discrete events can be compared quan¬ 
titatively with the variability-luminosity relationship us¬ 
ing Eq. (A6). But, it appears to us that discrete-event 
models do not satisfy the variability-luminosity relation¬ 
ship in a natural way. This difficulty possibly means that 
this kind of model does not apply to the active galactic 
nuclei, and, therefore, that Eq. (10) is not valid. Other 
source of variability (not considering the possibility that 
variations are extrinsic, which seems very unlikely) could 
be global effects, which do not have a simple general an¬ 
alytical formulation, in opposition to the local processes 
discussed here. 
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A. Appendix: Derivation of the variance of L{t) 


In our case, everything happens as if we have a super¬ 
imposition of an infinite (possibly continuous) number of 
shot noises given by the relationship: 

Lp ='^ ep{t - tip) (A3) 

ip 


(We neglect the constant term, because it does not appear 
in the variance). The p vector is the list of parameters that 
define the event; p is an element of the U(p) space, of 
arbitrary dimension, and follows a distribution P{p) with 
/■D(p) ~ 1- average rate of these shot noises are 

Np ■ dP{p). For each p, the variance of the shot noise is: 


dVar(Lp) = Np ■ dP{p) 


C + co 


ep{t) dt 


(A4) 


Therefore the variance of the total light curve is the sum 
(integration) of the variance of all shot noises: 


Var(L) = j dVar(Lp) = 

Jv(p) 


p p + ca 

/ Np - j ep{t) dt dP{p) 

IV(p) J-oo 


(A5) 


If one make the assumption that the rates Np are in¬ 
dependent of p, Eq. (A5) becomes: 


(•- 1-00 


Var(L) = A-/ / e|,(t) dt dP(p), (A6) 

-JT>{p) J-oo 


which gives Eq. (13). 


B. Appendix: Effects of the normalization and of 
the stretching factor on cr{L) 

Eq. (A5) can be developed further if we assign physi¬ 
cal meanings to some components of p. The normaliza¬ 
tion (i.e. the total energy liberated by an event) and the 
stretching factor are respectively the first and the second 
components of p. We can write: 


The variance of the generalized shot-noise process given 
in Eq. (10) can be calculated using Campbell’s theorem 
(e.g. Papoulis 1991), which states that if: 

s{t) =Y^ h{t-ti), (Al) 

i 

where the points U follow a Poisson distribution with a 
mean rate A, then the variance of s{t) is given by: 

f + co 

(A2) 


ep{t) =Pi - --ep ( — ) , (Bl) 

P2 ^\P2j 

where p = (pa,...). We impose ^p(t) dt = 1, so that 
we have Ep = ep(<) dt = pi, Vp. The reason why 

epi^) is multiplied byp^^ is that it makes ^p(t) dt 
be independent of p 2 (because f{k • x) dx = • 

f-^ /(cc) dx). Let us make the additional assumption 
that Pi and p 2 are independent of each other and of the 
other Pi’s. Let us write Vt oo), T>( ,e)> and Vip) the 


Var(s) = A • 


h'^{t) dt 
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spaces from which are drawn pi, p 2 , and p respectively. 
This makes Eq. (A6) become: 


Var(L) = N ■ 


f + CO 


1 


/ / Pi • — • 

'v(p) J-oo P2 ^ 


dt dP{p) = 


= N' 


f I f 

lv( ) Jt>( g) J-Dip) d-00 P 2 


V 


V" 


P\P2 


dt dP{p) dP{p2) dP{pi) = 


r+00 


= N-pI 


lv{ g) P 2 Jv(p) J — 00 
V 


P\P2 


dt dP{p) dP{p2) = 


r+00 


= N-pI 


Ivi e) ^2 Jv(p) J — 00 

V 


e-pit) dt dP{p) dP{p2) = 


= N-pIp^^ 


r+co 




This relationship is used to derived the scaling properties 
in Sect. 4.2. To obtain the results shown in this section, 
we introduce two random variables, Eq and with the 
relations: pi = E = k ■ Eq and p 2 = i = S ■ £q (we call K 
the multiple integral). The above equation becomes: 

YuriL) =N ■ (^S-^K (B3) 

With the same substitutions, the mean luminosity be¬ 
comes: 


L = N ■ {k-^ +C (B4) 

Eqs. (B4) and (B3) give directly Eqs. (15), (17), and (19). 
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